Tests for Gaussianity of the MAXIMA- 1 CMB Map 
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Gaussianity of the cosmological perturbations in the universe is one of the key predictions of standard 
inflation, but it is violated by other models of structure formation such as topological defects. We 
present the first test of the Gaussianity of the cosmic microwave background (CMB) on sub-degree 
angular scales, where deviations from Gaussianity are most likely to occur. We apply the methods 
of moments, cumulants, the Kolmogorov test, the test, and Minkowski functionals in eigen, real, 
Wiener-filtered and signal-whitened spaces, to the MAXIMA-1 CMB anisotropy data. We conclude 
that the data, which probe angular scales between 10 arcminutes and 5 degrees, are consistent with 
Gaussianity. 



A. Introduction: The anisotropy in the cosmic mi- 
crowave background (CMB) is arguably the cleanest cos- 
mic signal that preserves the intrinsic statistical proper- 
ties of cosmological perturbations The recent obser- 
vations from MAXIMA-1 § and BOOMERANG |] have 
unambiguously detected a sharp peak in the CMB power 
spectrum. This has favored inflation Q as the dominant 
mechanism for structure formation of the universe ^ as 
opposed to other candidates like topological defects [||. 
Another key prediction of standard inflation is that the 
distribution of cosmic perturbations are Gaussian, while 
other cosmological models such as isocurvature inflation 
(e.g. 0) and topological defects (e.g. ^) predict other- 
wise. Thus tests for the Gaussianity of CMB anisotropy 
data can discriminate between cosmological models. In 
addition, a Gaussian distribution is an important ingre- 
dient in the estimation algorithm of CMB power spec- 
tra |^,|l^, which has been used to produce the recent 
MAXIMA-1 and BOOMERANG results [|,|. Adding 
the fact that the estimation of cosmological parameters 
pT|-p^ relies on the estimated CMB power spectra, it 
is important to verify the Gaussian distribution of the 
CMB. 

Tests for the Gaussianity of CMB data have been 
carried out by numerous authors, mainly using data 
from the Differential Microwave Radiometer (DMR) on 
the COsmic Background Explorer (COBE) Q. Sev- 
eral statistics were used including moments, cumulants, 
Minkowski functionals (which include genus, e.g. ||l5| , |l6|| ), 
the three-point function (e.g. p^ ), bispectrum (e.g. 
|^,|9|), wavelet transform (e.g. ]20|]), etc. All these tests 
showed that the data were consistent with Gaussianity, 



except for a couple of results |jT^,^ that may have non- 
cosmological origins . 

In addition to the usual limitations from foreground 
contamination and instrumental noise, the 7° angular 
resolution of the DMR data is not ideal for tests of Gaus- 
sianity. Angular resolution is an issue because the size of 
the causal horizon at last scattering is about one degree. 
Thus in a sky patch of super-degree size, there are many 
uncorrelated perturbations due to causality. As a result, 
the central limit theorem guarantees that the pre-last- 
scattering anisotropy on super-degree scales will tend to 
be Gaussian. Even if the post-last-scattering anisotropy 
is strong and non-Gaussian, it can be obscured by the 
pre-last-scattering contribution, which we expect to be 
Gaussian. Therefore a CMB map with sub-degree reso- 
lution can provide much more powerful tests for Gaus- 
sianity. Park et al. |^ recently used a genus test on 
the QMAP and Saskatoon data, which have a resolution 
of about 1.5°, to show that they were consistent with 
Gaussianity. In this paper we report results from a series 
of Gaussianity tests on the MAXIMA-1 CMB map [||, 
which provides anisotropy information on angular scales 
between 10' and 5°. To optimize both the resolution and 
the signal-to-noise ratio for these tests, we shall use a 
map with 5972 square pixels of 8' each Q . Using these 
high-quality data we probe for the first time the Gaus- 
sianity of CMB anisotropy on sub-degree scales. 

A CMB sky is Gaussian if it is a realization of a pro- 
cess that is only specified by the two-point correlation 
function of a given cosmological model. However, be- 
cause we have only one realization of the CMB sky, it 
is not possible to exhaustively test the statistical prop- 
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erties of the process that generated it. Therefore, we 
employ a Frequentist approach, testing the null hypoth- 
esis of Gaussianity for the inflationary model that best 
fits the data. 

B. Karhunen-Loeve transform: We consider the 
Karhunen-Loeve (K-L) transform, sometimes called Prin- 
cipal Component Analysis or the signal-to-noise eigen- 
mode transformation [§^,D, because it enables us not 
only to transform the observed CMB map into uncor- 
rclated eigenmodes of known signal-to-noise ratios, and 
further to implement Gaussianity tests on the uncorre- 
lated modes. For the CMB, it is standard to model 
the data d = as a linear sum of uncorrelated signal 
s = Si and noise n = n^, with the correlation matrix 
C = (ddT) = S + N where S = (ss'^) and N = (nn^). 
In the noise-whitened space d*^^^ = N^^/^d, all the 
eigenvalues of the noise matrix, N^^/^NN^^/^ = I, 
are simply unity. Thus the eigenvalues e^^^^^ of the 
noise-whitened signal matrix N~^/^SN~^/^ represent the 
square of signal-to-noise ratios of each eigenmode. The 
coefficients b^^^ of the noise-whitened eigenmodes in a 
data set can be obtained by transforming d^^^ to the ba- 
sis which diagonalizes N~^/^CN^^/^. These coefficients 
are normally called the K-L coefficients. 

We compute S using the CMB power spectrum of 
the best estimated cosmological model (fif, = 0.105, 
ilc ~ 0.595, flA — 0.3, and h — 0.53) and including the 
effects introduced by the beam shape and the pixelization 
of the map [ p5[ . The matrix N is estimated from the tem- 
poral data p3[. The resulting eigenvalues in the noise- 
whitened space are shown in Fig. |l|, sorted in descending 
order. The dot-dashed line indicates that only the first 
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639 modes have signal-to-noise ratios [e'^'^^^] > 1. 
This number is well below our pixel number 5972. It is 
determined by the signal and noise levels and the observ- 
ing resolution (beam) of a data set, but is independent 
of the pixel number of the map when the pixel size is not 
larger than the observing resolution. 

A common technique to remove the information of 
non-cosmological interest in the map is Wiener filtering, 
jjWF _ sc-id. This is equivalent to weighting the eigen- 
modes with the ratios e'S'(w) /(e'S'CW) + i) (ggg pjg The 
sum of these ratios for all the eigenmodes is 837, again 
well below the total number 5972 of the eigenmodes. We 
shall employ this technique in section D to amplify the 
statistical significance of the CMB signal in our map. 

The K-L transform can also be used to test for 
Gaussianity. If the underlying map is Gaussian, 
then the eigenvalue-normalized K-L coefficients a'^^ = 
b(w) /(e'^(^) -1-1)^/2 should be a set of Gaussian variables 
with mean zero and variance one. We compute the a^^' 
of our data, and use the observed frequency distribu- 
tion to approximate its one-point probability distribution 
functions (PDF) p{i') , where v is the number of standard 
deviations from the mean, for both the entire 5972 and 
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FIG. 1. Noise-whitened eigenvalues of the MAXIMA-1 
CMB data (top), and the Wiener filter (bottom). 
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FIG. 2. PDF piv) of the entire 5972 (a) and the first 639 
(b) a'^' of the MAXIMA-1 data (solid), the Gaussian expec- 
tation (dashed), and its 95% CRG (shaded). 

the first 639 modes (i.e., those modes with signal domi- 
nating over noise). Both cases easily pass the and the 
Kolmogorov tests for Gaussianity at 95% confidence. 

In addition to the above tests, we also use Monte Carlo 
(MC) simulation to build a Gaussian reference frame for 
the Frequentist approach. We generate 100,000 Gaussian 
realizations, each of which is obtained by df^'^ — ClJ^gj, 
where gj is a Gaussian variable of mean zero and variance 
one. As a first application, we use the MC simulation to 
find the probability distribution of p(i^) at each for the 
entire 5972 and the first 639 a'^' . Fig. || shows that in 
both cases the real data lie well within the 95% confidence 
regions of Gaussianity (hereafter CRG). We also compute 
the moments and cumulants of the first 639 a*^^) up to 
tenth order, and they are all well within the 99% CRG. 
All these results support the conclusions not only that 
our map is consistent with Gaussianity, but also that our 
estimations of noise and CMB power spectrum (giving N 
and S respectively) are consistent with the data so as to 
provide the proper eigenmodes for the K-L transform. 
C. Minkowski functionals: The concept of Minkowski 
functionals is based on integral geometry. According to 
Hadwiger's Theorem /3-t-l Minkowski functionals are suf- 
ficient to measure the morphology of a /3-dimensional 
pattern. In the case of CMB we have (3 — 2 and thus 
need only three Minkowski Functionals. First we de- 
fine the excursion set Q in a given CMB field as the 
region in which the CMB amplitude di is larger than 
a threshold: Q = Q{v) — {di\{di — ^)/a > i/}, where 
/i = {d) and (j"^ — {d^) — /i'^. Then the surface densities 
Viiy) = Vi{Q)/A of Minkowski functionals Vi{Q) for a 
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FIG. 3. (a) Vertices of type I and II formed by inducing an 
offset in the pixel positions every other row. (b) An example 
of a pixelized CMB map with the artificial pixel offset. 



CMB patch of angular area A can be defined as 
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where dQ is the boundary of the region Q, dA and dl 
are the differential elements of Q and dQ respectively, 
and K is the geodesic curvature of dl. These Minkowski 
Functionals have different morphological meanings: Vo 
is the total area of Q, Vi is the total length of its 
boundary, and V2 is the number of isolated regions (hot 
spots) in Q minus the number of holes (cold spots). 
For an isotropic Gaussian field these functionals are 
characterized only by the field's variance and the 
variance of its gradient t = (|Vc?p)/2, i.e. Wo(g) — 
erfc(i//V2)/2, «i(g) = t1/2 exp(-i.V2)/8a, and W2(G) = 
rzyexp(— j/^/2)/(27r)^/^cr^. We also note that the com- 
monly used one-point PDF p(i^) and genus g, or the 
'Euler-Poincare characteristics', are simply related to the 
Minkowski Functionals. For the two-dimensional CMB, 
we have plv) = —dVo/Adv and g = ¥2 + Vb/27r. 

Because the MAXIMA- 1 map has an irregular bound- 
ary and all pixels have a square shape lying on a regular 
lattice, we need a numerical scheme to approximate its 
Minkowski Functionals. For a given threshold v and thus 
a given Q on the lattice, uq is simply the number of pixels 
in Q divided by the total number of pixels in the map. 
Here vi is determined by the length of dQ multiplied 
by 7r/4 to correct for the square-pixel lattice effect. To 
compute V2, we first induce offset in the pixel positions 
every other row. As Fig. || (a) shows, this gives two types 
(I and II) of vertices. For each type, the integration of 
K in equation (0) can be easily obtained by examining 
the status of the neighboring four pixels (basically the 
integration is determined by the change of the internal 
angle along dQ). Fig. ^ (b) shows an example, where the 
pixels in Q are shaded and the one with no CMB data 
is labeled with a cross. In this case the vertices ai, a2, 
03 and 04 have unambiguously {2K/T:)dl = —2, 2, —1 
and respectively, while a confusion in this calculation 
would occur at aiand ck2 if^ there was no offset every other 
row. We also notice that any bias possibly induced by 
the above algorithms should disappear when the number 
of pixels is large, as in our case. 

We compute the vo{iy), vi{iy), and ^2(1^) for both the 
entire map and the central 37% (2209 pixels, covering 
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FIG. 4. Minkowski functionals (solid) of the entire MAX- 
IMA-1 map (a) and the central part (b). The Gaussian ex- 
pectation obtained from the MC simulation (dashed), its 95% 
CRG (shaded), and the analytic Gaussian forms (dotted). 



about 6.3° X 6.3°), expected to have the lowest noise. We 
use pixel sizes of 8', 16', and 24', the last two obtained 
by averaging the neighboring pixels of the original map. 
The results of all six cases lie within the 95% CRG ob- 
tained from the MC simulation described earlier. Fig. ^ 
shows the results of the two cases with 8' pixels. We 
note that while the means of the MC simulation (dashed 
lines) are close to the analytical isotropic Gaussian forms 
(dotted lines) in (bl-b3), they deviate significantly from 
each other in (al-a3). This is due to the higher noise level 
near the edge of the map, contributing as an anisotropic 
component in the map (the RMS noise of the full map 
and the central part are about 162iiK and 57 respec- 
tively) . 

We also note that some results (solid lines) appear to 
have systematic departure from the Gaussian expecta- 
tions vf^'~^{i') (dashed lines). For example, (b3) shows 
the number of cold spots is less than the Gaussian ex- 
pectation. To determine whether these discrepancies are 
really systematic, we first define 



(^') 



dv', i = 0,l,2, (2) 



where a^'^{v') is the standard deviation of Vi{v') esti- 
mated from the MC simulation. Applying this to both 
the real data and the MC simulation, we obtain respec- 
tively the and its Gaussian-expected PDF p{Ii{v)). 
The results from all six cases lie well within the 95% 
CRG. Thus we know that the apparent systematic devia- 
tions from Gaussianity in Vi are statistically insignificant 
as measured by Ii{v). 

D. Wiener filtering and signal- whitening tech- 
niques: To increase the statistical significance of the 
CMB signal in the map, we now consider two filter- 
ing methods. One is the Wiener filtering addressed 
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FIG. 5. The Wiener-filtered (left) and the signal-whitened v v 

(right) MAXIMA-1 map, with ranges of fluctuations of FIG. 6. Similar to Fig. ^ but for the Wiener filtered (a) 

(-3.1,3.4)a and (-3.7, 4.3)cr respectively. and the signal-whitened (b) MAXIMA-1 central map. 



earlier, and the other is a new signal- whitening tech- 
nique = S^''^C~^d [Q, which not only removes the 
anisotropy on scales where the noise dominates (as in 
Wiener filtering) but also equalizes the anisotropy am- 
plitudes on scales where the CMB signal dominates. The 
signal- whitening procedure will reveal in the features 
of the non-Gaussian components whose contribution in 
the CMB anisotropy dominates the Gaussian one within 
at least a range of the accessible scales We apply 
these filtering methods to both the real map (Fig. ||) and 
the MC simulation, and then compute their Minkowski 
functionals. We find that for both the entire map and 
the central part (as previously) the Vi{i^) and /i(^) of the 
filtered maps are within the 95% CRG. Fig. ^ shows re- 
sults of the filtered central maps. We also verify for both 
of the filtered maps that none of the pixels with ampli- 
tude outside the ±2(t range coincides with the locations 
of any known radio or IRAS point sources Thus 
we have no statistically significant detection of localized 
non-Gaussianity in our data. 

E. Conclusion: We employ moments, cumulants, the 
Kolmogorov test, the test, and Minkowski functionals 
in eigen, real, Wiener-filtered and signal-whitened spaces 
to implement a total of 82 (not independent) hypothe- 
sis tests for Gaussianity (22 in Sec. B, 36 in Sec. C, 24 
in Sec. D), and show that the MAXIMA-1 CMB map is 
consistent with Gaussianity on angular scales between 10' 
and 5°. This gives us confidence in the Gaussian distri- 
bution used in our estimation of the CMB angular power 
spectrum 1^ , and consequently in the estimation of cos- 
mological parameters |ll|Jl3|1 . Although our results are 
consistent with standard inflation and the CMB power 
spectrum of recent observations has favored inflation as 
the dominant mechanism for structure formation , the 
existence of topological defects is still possible. More so- 
phisticated methods or data with even higher resolution 
and signal to noise ratio will be required to fully explore 



this possibility. 
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